Suppression of carrier induced ferromagnetism by composition and spin fluctuations 

in diluted magnetic semiconductors 
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We suggest an approach to account for spatial (composition) and thermal fluctuations in "disor- 
dered" magnetic models (e.g. Heisenberg, Ising) with given spatial dependence of magnetic spin- 
spin interaction. Our approach is based on introduction of fluctuating molecular field (rather than 
mean field) acting between the spins. The distribution function of the above field is derived self- 
consistently. In general case this function is not Gaussian, latter asymptotics occurs only at suffi- 
ciently large spins (magnetic ions) concentrations rij. Our approach permits to derive the equation 
for a critical temperature T c of ferromagnetic phase transition with respect to the above fluctua- 
tions. We apply our theory to the analysis of influence of composition fluctuations on T c in diluted 
magnetic semiconductors (DMS) with RKKY indirect spin-spin interaction. Two new results are 
obtained, (i) It is shown that composition fluctuations destroy ferromagnetic ordering in DMS if the 
carrier concentration n e exceeds the critical value n c « O.lrii.(ii) In the case n e < 0.5n c and mag- 
netic ion spin S = 5/2 the composition fluctuations shift the critical temperature T c (as compared 
to its mean field value) by a factor about I — 5n e /ni. 



Ferromagnetic (FM) ordering in p-doped diluted mag- 
netic semiconductors (DMS) attracts much attention of 
both experimentalists 0- || and theorists [||- Com- 
mon approach to this problem is a mean field approxima- 
tion (MFA) , corresponding to some mean overall (homo- 
geneous) carrier and magnetic ion magnetization. Once 
appeared, the exchange fields of carriers and magnetic 
ions lead to mutual spin splitting that can be stable at 
low enough temperatures T < T c The problem of 

"high temperature ferromagnetism" in DMS is now very 
important from the point of view of technical applica- 
tions. So, the question how to increase the critical tem- 
perature T c has been widely discussed in recent publica- 
tions (see H and references therein). One way to control 
T c is to vary the concentrations of magnetic ions and/or 
carriers. In the case of degenerate carriers, the critical 
temperature T c of FM phase transition in the bulk DMS 



1 was pointed out earlier ||, but quantitative contribu- 
tion of disorder in the magnetic ion subsystem has not 
been considered till present. Moreover, the problem of 
the existence or non-existence of a finite critical temper- 
ature T c > at arbitrary relation between rij and n e has 
not been resolved also. 

Here we develop a theory of magnetic ordering of dis- 
ordered magnetic system with given spin-spin interaction 
in terms of fluctuating local molecular field approxima- 
tion (FFA). The method we propose here is applied to 
calculations of T c caused by RKKY interaction. 

The Hamiltonian of the problem reads 



(i) 



3<3 



is proportional to n,n, 



1/3 



where rii and n e are the 



magnetic ion and carrier concentrations respectively. 

Beyond the MFA, we should consider the Friedel os- 
cillations of carrier spin polarization Q occurring at the 
scale of inverse Fermi wave vector 1/kp. In this case, the 
fluctuations of inter-magnetic ion distance r<y should be 
taken into account, if the nearest neighbor mean value r 
is about 1/fcir. Really, in the case of kpf 1 the num- 
bers of magnetic ion pairs with FM and AFM RKKY 
interaction is comparable, which can lead to suppression 
or complete destruction of long range FM order. 

Note, that the limitation of MFA by inequality kpf << 



where magnetic field Hq and interaction J(ty,-/) is mea- 
sured in energy units (i.e. g/i = 1, fi is Bohr magneton). 
The Hamiltonian ([!]) incorporates two "sources of ran- 
domness". First, (the spatial disorder) is that spin can 
be randomly present or absent in the specific j-th cite of 
a host semiconductor. Second, (the thermal disorder) is 
a random projection of a spin (if any) in j'-th cite. 

Note, that our general formalism is valid for any form 
of Jifjji) so that its specific form will be chosen later. 

In a mean field approximation, the Hamiltonian (|l|) 
reduces to the sum of Zccman energies 



H = Y,(Hi + H )Si 



(2) 



1 



in the external magnetic field Hq and local molecular field 

where \ is a thermal average at a site j. Next step of 

the MFA is to substitute the molecular field @) by the 
mean field 



»iJ J 



(4) 



The bar means the averaging over spatial disorder , so 
the H m f is homogeneous over crystal volume; I^Sj is an 
average spin in this mean field, 



S) = 



Tr{Se-^f+ a o)S } 



(5) 



T re -0(H m j+H o )S 

(3 = 1/T. The system of Eqs (g) and (§) gives well- 
known solution for mean fields H m f and magnetization 

M = ng(j,(§y 

This work draws attention that the spatial and ther- 
mal fluctuations can be taken into consideration by in- 
troduction of random field rather than mean field. In 
our approach, we consider every spin Sj as a source of 
fluctuating (random) field 



Hf EE J(fi 



(6) 



affecting other spins at the sites ?i- In other words, ev- 
ery spin in our approach is subjected to some random 
(rather then mean) field, created by the rest of spin en- 
semble. So, all thermodynamic properties of the system 
will be determined by the distribution function f(Hf) of 
the random field Hf. Namely, any spin dependent macro- 
scopic quantity (like magnetization) << A » reads 



<< A >>= 



< A > S} f(H f )dH f , 



(7) 



where < A >g is single particle thermal average with 
effective Hamiltonian (Q), where Hi is substituted by Hf. 
The distribution function f(Hf) is defined as 



f(H f ) 




rj)Sj-H 



(8) 



Our FFA approach is based on micro-canonical statis- 
tical theory of magnetic resonance line shape ||. Lat- 
ter theory assumes the additivity of local molecular field 
contributions as well as the non-correlative spatial distri- 
butions of spins (magnetic ions). 

Latter assumptions with respect to spectral represen- 
tation of <5 function permit to transform the Eq. (^]) to 
the self-consistent integral equation for f{Hf) = f{H). 
Introducing the probability n,i(r)d r for small volume d 3 r 
to be occupied by a particle, we obtain 



f{H) = J exp \it{H - H Q ) 
where according to definition (0) 

<< exp[-iJ(r)5f] - 1 >>= 



exp ( / << exp 



iJ(f)Sf - 1 >> rii(r)d 3 r 



d 3 T 



Tr 



exp 



(-iJ(r)§Tj exp(-/3(i? + H )S) 



Tr 



exp 



f(H)d 3 H - 1 



(9) 



(10) 



Eqs (|9|) and (|10|) represent the integral equation for 
distribution function f(H). This equation can be solved 
iteratively. 

However, in many cases it is possible to avoid the so- 
lution of the integral equation since in these cases it is 
exactly reducible to the set of transcendental equations 
for macroscopic quantities like << S » (magnetiza- 
tion),..., << S n >>, n > 1 of the system. 

This approach can be most easily demonstrated for the 
Ising case in Hamiltonian (Q). Although this case is less 
general than Heisenberg one, it is perfectly valid either 
for the case of hole induced ferromagnetism in quantum 
wells or in the case of host semiconductors with uniaxial 



symmetry. 

If the magnetic fields are directed along OZ axis, the 
scalar product reduces to St = Szt = tut, m — —S, ...,S 
and (S) takes the form 



f(H) = j- I e tHT e g ^dT- 



(11) 



Q{t) = ((f "iff (e- u ^ mT - l) d 3 ^. (12) 
The calculation of average in (ll^) yields 
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e -iJ(7)mT - A\ = K [cos( J(r)mT) - 1] + ib m sin( J(f)mr)} . 



m=l/2 



a m = f A m (/3H)f(H)dH, b m = J B m ((3H)f(H)dH. 

A m (/3H) = o C0S Mm/3ff)sinh(/3ff/2) ^ B , m) = ^(mpH) sinh(/3ff/2) 



sinh((S + 1/2)0H) 



sinh((S + 1/2)PH) 



(13) 



In this case Eq.([l2]) assumes the form 

<7(t) = £({a m ,6 m } ,r) = 
s 

= {a m fi(mr) + ib m T 2 {mr)} ; (14) 

m=l/2 

^i(mr) = / rii(r) [cos( J(f)mr) — l]d 3 r, 



fjjmr) = / ni(r) sin (J (f)mT)d 3 r. 



(15) 



Let us return to Eq.(^) for f(H). Multiplying it by 
A m ((3H) and B m (j3H) and integrating over H, we obtain 
the closed system of nonlinear equations for the numbers 
a m and b m . The explicit form of this system reads 



b m = 



A m {T,f3) 



f * 








| cos 


^ b ml T 2 {m'T) dr; 


^ro'=l/2 , 




\m'=l/2 / 






/ S \ 


a m >T\{mT) 


> sin | 


^ b m ,Ti(m!r) dr, 


,m'=l/2 J 




\m'=l/2 J 


(ffr)dff; B m (T,j9) = 


4i 


fOO 

B m (0H)sm(Hr)dH. 

3 



(16) 
(17) 



To obtain (|16j), we substite of (15j) into Eq.(y 



pears. In terms of Eqs (|16[) this means that we can 



Solutions of the Eqs (|16| ) allow to obtain magnetiza- put b m — > for T > T c . In this case, the parameters 
tion of the system .M = gfxMi, where Mi = — ((Sz)) = 



Em=i/2 TO ^m as weu as other spin averages 



M k ee = {-if £ m h 

m=l/2 



a m ; for even k, 
bm ; for odd /c. 



(18) 



The critical temperature T c is usually defined as a tem- 



a m = a m (0) should be found from the equations 
a m = J A m (T,{3)exp < a m ,T\{m!T) > dr; 



= 1/2 



(19) 



that stem from ( |16| ) for fe m = 0. 

Parameters 6 m can be obtained by linearization of the 



perature, where nonzero infinitesimal magnetization ap- second equation (hq) 



bra — ^ ^ Km.m'brn 1 \ K-w 
m'=l/2 



/ B m (T, /5)J r 2 (m , T) exp 

J -°° [m"=l/2 J 



(20) 
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Thus, near T c parameters b m obey the system of homo- 
geneous linear equations. Well-known condition of exis- 
tence of non-trivial solution of this system gives following 
equation for T c : 

det|K-I|=0, (21) 

where I = 5 m m i is the identity matrix, m, ml = 
1/2,. ..,5. 

The Eqs |f|) and are the main result of the- 
oretical background we develop. We start the analy- 
sis of these equations from the case of S = 1/2. The 
system (^0|) reduces to single equation with ai/ 2 = 
A l/2 = 1;A/ 2 (t,/3) = 5(t) and B 1/2 = tanh/?ff/2; 

Subsequent calculations require the definition of spa- 
tial dependence of J(f). Usually in the problems of car- 
rier -induced ferromagnetism in DMS, the RKKY in- 
teraction Q is considered as an effective spin-spin ex- 
change interaction that results in FM ordering. Re- 
cently, the complex valence band structure of III-V and 
II- VI semiconductors was discussed as a reason for mag- 
netic anisotropy and some enhancement of carrier in- 
duced magnetic interaction JlOj] , ||. To clarify the role 
of fluctuations, here we neglect the effects of complex va- 
lence band structure, which can be also incorporated in 
our theory. 

In the case of simple one band carrier structure, the 
RKKY interaction reads 



3 MF 



J{r) = J a - 



Jo 



x = 2kpr 



47T 3 



(22) 
(23) 



where J C i is a carrier-ion exchange constant, fio is a unit 
cell volume, m* is the density of states effective mass. 
The critical temperature is determined now by the equa- 
tion that follows from the Eq. (]2C|) after the substitution 



of Eq. (22) to integrals (15 



1 



(?) 

2 y^sinh^) 



exp 



ip! (£) d£. (24) 



Here we introduce the dimensionless parameters 8 — 
ttT c /2Jq and v = rii/6ir 2 n e . The functions ipi (£) and 
¥>2 (?) are related to integrals ( [l5| ) . In the case of homo- 
geneous magnetic ions distribution, rii(f) — m = const, 



<pi (0 

(0 



1 



.XCOSX — SIM 



x cos x 



OCT* 

sinx 



x dx 



x dx. 



(25) 
(26) 



The MFA results (for spin S = 1/2) can be recovered 
from Eqs ( p5| ) and ( p6| ) by their expansion up to the linear 
terms ipi (£) — » 0; (f2 (?) — > £, that leads to the equation 
for critical temperature in the form 



V7T 2 /4. 



Gaussian asymptotics of distribution function can be 
used to analyze the influence of fluctuations on T c near 
MFA. Latter asymptotics corresponds to the next term 
of expansion of the Eq. (|2l! 



(27) 



<Pi (?) - < 2 /30, Vi (0 - i- 



The numerical solution of the Eq. (|I|) with respect 
to ( |27| ) is presented in the Fig.l. Note the qualitative 
difference in behavior of 9 as a function of n e /rii for 
MFA and Gaussian distribution function. While MFA 
predicts the FM ordering at any relation between the 
electron and magnetic ion concentrations, the fluctua- 
tions suppress ferromagnetism and make it impossible if 
the concentration ratio n e /rii = 1/0)Tt 2 v exceeds critical 
value {n e /rii) c = 1/Qtt 2 v c . We find 



(^)Gauss = vr/15 w 0.2094 



(28) 



as a solution of the Eq. (|24|) in the limit 9 — > that corre- 
sponds to fluctuation restriction for carrier concentration 

(Kc)Gauss = 5/27T 3 « 0.08^. 




FIG. 1. Dimensionless phase transition temperature 
T c I Joo for spin 1/2 versus ratio v a / Vi. Here 



(6tt 2 ) 4 / 3 J^mTj 2 / 3 



n >± - — , Ui >e = n.i, e flo, fio is unit cell volume. 

Here we put Vi = 1. 1- non-Gaussian distribution function, 
2- Gaussian distribution function, 3 - MFA. 
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We can also take into account the deviations of molec- 
ular field fluctuations from their Gaussian asymptotics 
by numerical solution of the Eq. (|^). The results are 
also reported in the Fig.l. The critical concentration for 
this case turns out to be 



0.1366, 



0.12m 



The visible deviation of T c (|24|) from that for Gaus- 
sian asymptotics is revealed at relatively narrow interval 
n e /rii near the critical value (n e /rii) c . It can be shown 
that this property holds for any other magnetic ion spin 
5 > 1/2. This permits to assert that the molecular field 
fluctuations can be considered with sufficient accuracy 
with Gaussian asymptotics for distribution function. 

Gaussian asymptotics for 5 > 1/2 corresponds to the 
expansion of Eq. ( |l5| ) in powers of t up to the second 
order (it can be shown that this is valid for sufficiently 
high spin concentration): 



with 



t- I \ 2 2 2 

>i(mr) = —a m r , 
.7-2 (tot) = —AHuit 



rii{r)J 2 {r)cPr, 



AH = I m(r)J(r)d a r. 
'v 



(30) 

(31) 
(32) 



Approximation (|30|) permits to find the explicit form for 
Gaussian asymptotics of distribution function ( pi] ) . This 
function depends on only two parameters ([l8]) Mi and 
M 2 : 



fc{H) - fM u M 2 (H) 

1 



2ctVttM2 



exp 



(H - AH Mi) 
4ct 2 M 2 



(33) 



Thus, instead the (25+l)/2 equations (16) the Gaussian 
asymptotics " generates" only two of them 



Mi = 



(Sz) Hf fM u M 2 (Hf)dHf, 

/oo 
(S%) f Ml ,M 2 (H f )dH f . 
-oo 



(34) 
(35) 



Note that the equations ([34j) and (^) resemble replica- 
symmetric solution found for Ising (spin 1/2) spin glass 
(see Q). The reason we do not use replica formalism 
here is that the explicit form of J(r) is important for 
DMS. The MFA result immediately follows from the 
Eq. (34) in the limit of zero fluctuations, a — > 0, 
f Mu M 2 (H) ^ 5(H - AHMi): 



T, 



MF 



-5(5 + 1) / ni (r)J(r)d 3 r. (36) 
3 Jv 



The equation for T c can be obtained from Eq's ( |34| ) 
and (^||) as Mi — > 0. Then Eq.(Q) can be transformed 
to following equation for critical temperature 



AH 



2a^JM 2 



dx, 



(37) 



(29) where mg (x) — SBs(Sx) is non-normalized Brillouin 



function, m' s (x) is its derivative. In the case of small 
fluctuations, 2ay/M^ «T, Eq. @ is reduced to MFA 
result @. 

We can find analytically the first fluctuation correction 
to T^ F due to Gaussian fluctuations of molecular field. 
The first term of expansion of (B7f in a / AH yields 



rpMF 



1 



1 



1 



AH 2 



(38) 



25(5+1), 

The physical meaning of Eq.(|38|) is the lowering of T c 
by the fluctuations. This means that due to fluctuations, 
FM ordering occurs at lower temperature or that the fluc- 
tuations suppress FM order at a given temperature. 

To apply this result to the case under consideration, 
the parameters (pl|) and (p2\) have been evaluated for 



(39) 



RKKY interaction (p2[) and constant spins concentration 
rij(r) = rii = const: 

Jo rij 

127T Tie. 

Jo 



AH = 



a = 



This gives 



rpMF 



= 1 - 



127+ 

~25~ 



1 



1 



25(5+1) 



Tie 

m 



RKKY 

For the important case of the magnetic ions Mn 2+ 
5 = 5/2 the Eq. © yields 1 - 5.0(n e /n l ). 



(40) 
with 




v /v. 
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FIG. 2. Dimensionless phase transition temperature 
T c /Jqq for spin 5/2 versus ratio v e /vi (also Vi — 1). The 
notations are similar to those for spin 1/2. 1- Gaussian dis- 
tribution function, 2- first fluctuation correction to MFA, 3 - 
MFA. 

In Fig. 2, we reproduce the electron concentration de- 
pendence of T c calculated in MFA, with Gaussian dis- 
tribution function (|37j) and with approximation ( ^0| ) for 
S = 5/2. The dependence T c — T c (n e /rii) presented 
in the figure reveals also the presence of critical ratio 
(n e /rii) c . The equation for the {n e /rii) c can be found as 
T -+ in Eq.@. Equation we obtain {{AH) 2 = ira 2 ) 
is independent of spin, so account for the Eqs ( |39| ) re- 
produces the result we discussed for spin S = 1/2: 
v c = tt/15, or {n e / ni ) c = 5/2tt 3 ~ 0.08. 

In conclusion, here we present a theory of fluctuated 
molecular field for the spin systems with given depen- 
dence J{r). We derive the equations for order parame- 
ters and FM phase transition temperature for Ising model 
with indirect spin-spin RKKY interaction in DMS. Our 
analysis show, that fluctuations of transverse spin com- 
ponents (which are present in Heiscnberg model) do not 
reveal new qualitative properties in the problem of FM 
ordering in DMS. 

Naturally, the RKKY interaction is not unique spin- 
spin interaction in DMS. There is additional background 
interaction independent of free carrier concentration. In 
II- VI DMS this AFM exchange interaction is dominant 
for close magnetic pairs. An intuitive approach to this 
problem was proposed in Ref fujl , where MFA was ap- 
plied to ensemble of magnetic ion spins with excluded 
AFM exchange pairs. Our formalism can naturally in- 
corporate this interaction by adding correspondent terms 
to Hamiltonian ([l]). 

Actually, the parameter of theory that present the fluc- 
tuations of RKKY interaction in bulk DMS is a cube 
of ratio of Fermi wave vector length to mean inter-ion 
separation. We can expect that similar parameter will 
be "responsible" for fluctuations in 2D systems. This 
mean that critical temperature would depend on carrier 
concentration in accordance with MFA prediction in 2D 
case. 

The detailed analysis of all aforementioned factors will 
be presented elsewhere. 



ACKNOWLEDGMENTS 



This work was partially supported by the Grant No. 
09244106 of the Ministry of Education and Science of 
Japan. 



[1] T. Story, R. R. Galazka, R. B. Frankel and P. A. 

Wolff , Phys. Rev. Lett. , 56, 777 (1986); M. Z. 

Cieplak, A. Sienkiewicz, Z. Wilamowski and T. Story, 

Acta Phys. Pol. A 69, 1011 (1986); T. Story, G. 

Karczewski, L. Swierkowski and R. R. Galazka, Phys. 

Rev. B 42, 10477 (1990); T.Story, P.J.T.Eggenkamp, 

C.H.W.Swuste, H.J.M.Swagten, W.J.M. de Jonge and 

A.Szczerbakow, Phys. Rev. B 47, 227 (1993). 
[2] H.Ohno, H.Munekata, T. Penney, S. von Molnar and 

L.L. Chang, Phys. Rev. Lett. 68, 2664 (1992). H.Ohno, 

A. Shen, F. Matsukura, A. Oiwa, A. Endo, S. Kat- 

sumoto, and Y. Iye, Appl. Phys. Lett. 69, 363 (1996). 

F. Matsukura, H.Ohno, A. Shen and Y.Sugawara, Phys. 

Rev. B 57 R2037 (1998). 
[3] E.A.Pashitskij and S.M.Ryabchenko, Fiz. Tverd. Tela, 

21 545 (1979) [Sov. Phys. Solid State, 21, 322 (1979)]. 
[4] J. Konig, Hsiu-Hau Lin, and A. H. MacDonald, Phys. 

Rev. Lett. 84, 5628 (2000). L. Brey and F. Guinea, Phys. 

Rev. Lett. 85, 2384 (2000). 
[5] C. Zener, Phys. Rev. 81, 440 (1950); 83, 299 (1951). 
[6] T.Dietl, H.Ohno, F. Matsukura Phys. Rev. B 63, 195205 

(2001). 

[7] Ch. Kittel Quantum Theory of Solids. , John Wiley & 
Sons, New York, 1987. 

[8] A.M.Stoneham, Rev. Mod. Phys. 41, 82, (1969). 
V.A.Stephanovich Ferroelectrics, 192, 29 (1997). 

[9] Corenblit and Shender, Usp. Fiz. Nauk, 157, 267, (1989). 
[Sov. Phys. Uspekhi, 32, 139, (1989)]. 
[10] M. Abolfath, T. Jungwirth, J. Brum, and A. H. Mac- 
Donald, Phys. Rev. B 63, 054418 (2001). 
[11] D.Ferrand, J. Cibert, A. Waisiela, C. Bourgognon, S. 
Tatarenko, G. Fishman, T. Andrearczyk, J. Jaroszynski, 
S. Kolesnik , T. Dietl, B. Barbara, and D.Dufeu, Phys. 
Rev. B 57 085201, (2001). 



6 



